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ABSTRACT 

We present a sample of 20 runaway M dwarf candidates (RdMs) within 1 kpc of the Sun whose 
Galactocentric velocities exceed 400 km s“^. The candidates were selected from the SDSS DR7 


M Dwarf Catalog of West et al. (20111. Our RdMs have SDSS+USNO-B proper motions that are 
consistent with those recorded in the FFMXL, LSPM, and combined IL7S'i?+SDSS+2MASS catalogs. 
Sixteen RdMs are classified as dwarfs, while the remaining four RdMs are subdwarfs. W e model the 
Galactic potential using a bulge-disk-halo profile ( [Kenyon et al.|20^ Brown et al.|2014 ). Our fastest 
RdM, with Galactocentric velocity 658.5 ± 236.9 km s“^, is a possible hypervelocity candidate, as it 
is unbound in 77% of our simulations. About half of our RdMs have kinematics that are consistent 
with ejection from the Galactic center. Seven of our RdMs have kinematics consistent with an ejection 
scenario from M31 or M32 to within 2a, although our distance-limited survey makes such a realization 
unlikely. No more than four of our RdMs may have originated from the Leo stream. We propose that 
to within measurement errors, most of our bound RdMs are likely disk runaways or halo objects, and 
may have been accelerated through a series of multi-body interactions within the Galactic disk or 
possibly supernovae explosions. 

Keywords: brown dwarfs - proper motions - stars: kinematics and dynamics - stars: late-type - stars: 
low-mass 


1. INTRODUCTION 

Recent studies have shown that not all stars within 


the Milky Way are bound to the Galaxy (Brown et 


al.| 120051 Ede 

mann 

et al. 

120051 Brown et al. 

jljOObalbl 

2007a|bl 12010 

120121 

:^ol4| 

Falladino et ai.||201^ 

fal IZhengl 

et ai.||20i4p. 4 

’auris ( 

^61^1 

defines genuine hypervelocity 


et al.j (|2005|) thus proposed a solution that this star is 
a blue straggler, a blue main-sequence star (such as B 
type) that persi sts in clusters beyond th e turn-off point 
of the cluster. [Edelmann et al. (20051 suggested that 


tial of our Galaxy.” Stellar escape velocity is a function 
of Galactic position. Th e escape velocity th reshold is 
^400 km s“^ in the halo (Kenyon et al. 2008), whereas 


the current measured escape vel ocity tmeshold near the 


Sun is about 550 ± 50 km s ^ (Smith et al. 


2007 


Piffl 


such stars may come into existence through the merger 
of two lower-mass stars. A related merger-type event 
between two helium white dwarfs has been proposed to 
explain the presence of the hyperveloc ity, subluminous 
0-type star US 708 ( Hirsch et al.||2005 l, in the Galactic 
halo, with a reported GC velocity of 708 ± 15 km s“^. A 
more-recent study proposed that the energy of a su per- 
nova propelled US 708 outward (Geier et al. 2015) and 


et al.|[20l4 ). Hypervelocity stars are important, because reported an updated G G velocity of 1157 ± 53 km s 


they, unlike the vast majority of stars in the Galaxy, are 
not bound. One or more mechanisms are thus required to 
explain their hypervelocity origins. Star s may reach hy¬ 
pervelocity speeds, as predicted by Hills (1988), through 
a gravitational interaction between a binary system and 
the central massive black hole (CMBH), in which one of 
the stars is ejected from the Galaxy while the other is 
captured into orbit around the CMBH . The hyperveloc- 


Palladino et al. (2014a) presented a sample of 20 G and 
K hypervefocity candidates whose kinematics were in¬ 
consistent with a GMBH ejection scenario. With the 
aid of impro ved proper motion measurements, jZiegerei] 


ity sta rs studied by Brown et al. (2014) and Zheng et al. 
(2014) have kinematics that are consistent with ejection 
radially outward from the Galactic center. 

Alternative mechanisms may be responsible f or ob¬ 


served high Galac tocentric (GG) stellar velocities. Edel¬ 


mann et al. (2005) studied the main-sequence B-type star 
HE 0437-5439, whose velocity is at least 563 km s“^ 


Edelmann et al. (2005) showed that if HE 0437-5439 has 
a GC origin, then its flight time from the center to its 
current location would be ~100 Myr, about four times 
longer than its age on the main sequence. jEdelmaim] 

^ Department of Physics and Astronomy, University of Maine, 
Orono, ME 04469 

^ Department of Astronomy, Boston University, Boston, MA 
02215 


et al. (2015) dismissed a CMBH ejection scenario for 14 
ot these stars and concluded that they are likely disk 
runaways. 

While O and B stars may remain on the main sequence 
for ~10® yr, the lowest-mass M dwarfs may rernain on 
the main sequence for ~10^^ yr (jLaughlin et al. 1997). 
As M dwarfs (dMs) comprise ^707o ot the stars in the 
Milky Way, they are the dominant stellar constituent 
( Bochanski et al.|[2010| . Thus, depending on the mecha- 
nism that accelerates these stars to hypervelocities, one 
might expect to find more hypervelocity M dwarfs than 
hypervelocity O or B stars. In the last few years, several 
new studies of hypervelocit y candida t es late r than B-type 


sta rs have been published. Li et al. (2012), Palladino et 


al. (2014a) and Zhong et al.'f(2014) present hypervelocity 


candidates of stars as late as K type. However, very few 
M type hyperveloci ty stars have been observed to date. 


Vickers et al. (2015) present 18 hypervelo city candidates 
ranging in types from F to M. In addition, jSavcheva et al.j 
















































































































2 


Favia, West, & Theissen 


(2014) includes 14 new hypervelocity M subdw arf candi¬ 
dates with GC velocities exceeding 525 km s“^. [Theissen 


& West (2014) also present a sample of hypervelocity M 
dwart candidates. 

The primary challenge when observing M dwarfs is 
their intrinsic faintness, which presents a major limi¬ 
tation for detecting and studying hypervelocity stars. 
Along the main sequence, B-type stars are ~ 10^ — 10"^ 


times more luminous than M dwarfs (Habets & Heintze 


1981). The low luminosity of M dwarfs also makes spec¬ 
tral and proper motion analysis of M dwarfs particularly 
challenging. The discovery of runaway M dwarfs would 
thus provide a significant contribution to the population 
of identified runaway stars and opens up the possibil¬ 
ity of finding hypervelocity M dwarfs. A larger sample 
of M dwarfs is needed to identify more potential hyper¬ 
veloci t y cand idates. The spectroscopic catalog of [West 


et al. (2011 hereafter Wll) consists of 70,841 visually 
inspect ed A1 dwarfs fro m the Sloan Digital Sky Survey 


(SPSS , Yor k et al. 2000) Data Release 7 (DR7, Abazajian 
et al. 2009[) . The kinematic data from the Wll catalog 


thus present an ideal starting point for identifying new 
hypervelocity candidates. 

In this Paper, we present a sample of runaway M dwarf 
candidates (RdMs) from the Wll catalog. We discuss 
the process by which we create a high-quality sample 
of RdM candidates in Section We present the stellar 
properties of our RdMs and determine the probabilities 
for each of our RdMs to be bound or unbound in Section 
We discuss selection effects and possible acceleration 
mechanisms for our RdMs in Section |4| We summarize 
our Paper in Section]^ 

2. DATA 
2.1. Background 

The stars in Wll we re inspected using the Hammer 
spectral typing facility ( Covey et al.||2007 ) to determine 
spectral types and ensure that extragalactic interlopers 
were excluded. Stars were included only if they had 
S/N > 3 at ~8300 A. Proper motions of the stars in Wll 
were ext racted from the SDSS-I-U SNO-B proper motion 
catalog (Munn et al. 2004 2008). Distances were cal¬ 
culated using the photometric parallax relation between 
Mr and the extinction-c orrected r — z color of the star 
(Bochanski et al. 2010). Heliocentric radial velocities 


(RVs) were determined by cross-correlating each spec¬ 
trum with a template for t hat subtype, yielding t ypical 
precisions of 7-10 km s“^ (Bochanski et al. 2007 here¬ 
after B07b). The M dwarfs in the Wll catalog were 
categorized based on whether they passed a set of pro¬ 
cessing Hags, described below. 

Proper motions are essential in measuring the tangen¬ 
tial velocities (TVs) of M dwarfs. The hypervelocity O 
and B-type stars previously studied are mostly RV val¬ 
idated, whereas nearby RdMs can have high RVs, high 
proper motions, or both. Proper motions for stars in the 
W ll catalog were pul l ed fro m the proper motion catalog 
of [Munn et al.j (j2004[ |2008|). We ad opted a set o f clean 
proper motion flags established by Kilic et al. (2006), 
which enabled large, statistically complete sainples of 
white dwarfs from SDSS. Several of their processing flags 
include: match, the number of objects in USNO-B which 
match an SDSS object in a 1” radius; sigRA & sigDec, 


the residuals in the linear fit (based on all epochs at 
which the object was observed) of right ascension and 
declination proper motions; nFit, the number of detec¬ 
tions between SDSS and USNO-B used in the fit, with a 
maximum of 6; and dist22, the angle (in arcsec) to the 

nearest neighbor in SDSS with g < 22. _ 

The distribution of proper motion errors in jMunn et al.j 
(2004), who required match = 1, sigRA & sigDec < 350, 
and nFit > 4, can b e assessed by exam ining the proper 
motions of quasars. Munn et al. (2004) state that “ob¬ 
jects detected on fewer than tour plates in USNO-B suffer 
fr om a large contam ination by false matches.” According 


to Dong et al. (2011), contamination in the sample can be 
reduced by requiring match = 1, sigRA & sigDec < 525 
(the condition for the residuals is relaxed), nFit = 6, 
and dist22 > 7. The condition nFit = 6, according to 
Dong et al. (2011), “makes the cleanness of our samples 
much better than that of samples defined us ing the stan¬ 
dard criteria defined by Munn et al. (2004), even when 
we loosen the requirement on rms httmg residuals a lit¬ 
tle bit.” For comparison, the processing flags adopted by 
Wll include: GOODPHOT, which equals 1 if the distances 
are defined by the photometry and have measured mag¬ 
nitudes, 0 otherwise; GOODPM, which equals 1 if the proper 
motions satisfy the conditions match = 1, dist22 > 7, 
sigRA & sigDec < 1000, and (nF it = 6 or (nFit = 5 and 
(O < 2 or J < 2)))|Munn et aL] (|2004|), 0 otherwise (see 
footnote 13 in West et al.||2011 ); and WDM, which equals 
1 for a possible white dwart-M dwarf pair, 0 otherwise. 

Proper motions in the SDSS-FUSNO-B catalog com¬ 
bine relative proper motions from observations taken on 
USNO-B plates with SDSS CCD photometry. Stellar 
proper motions are also recorded in other catalogs, such 
as the Positions a nd Proper Motions catalog (PPMXL 
Roeser et al. IMot. the Lepine Shy a Proper Motion cat- 
alog (LSPM, Lepme V Shara|2005 ), which was inspected 
by eye, and the Wide-fie la Infrared Survey Explorer 
( WISE, [Wright et al.||2010|) , which is tied to the Inter¬ 
national LlbbrdmateRel^nce System (IC RS) through 
the Two -M icron All-Sky Surv ey (2MASS, iSkrutskie et 

:cf 


al. 2006). Dong et al. (2011) investigated the proper 


motion error distribution of the SDSS-I-USNO-B catalog 
using quasars, as they have extragalactic distances and 
so are expected to exhibit no observable proper m otion. 
Using the 1 0* ^ release of the SDS S quasar catalog (Paris 


et al. 2014), Dong et al. (2011) found that the proper 
motion errors of quasars are well represented by a dis¬ 
tribution consisting of a Gaussian core with extended 
wings. 


The catalog by Theissen et al. (submitted) established 
updated proper motions for > 10® KML stars from a 
combination of SDSS, WISE, and 2MASS baselines. The 
choice to adopt a baseline depended on whether or not 
each source passed a set of processing flags. Within their 
catalog, processing flags from among the SDSS, WISE, 
and 2MASS catalogs were applied to the SDSS quasar 
catalog to define a set of QSOs. These were used to 
compute angular distances between catalog baselines and 
provide updated proper motion errors for KML stars, 
depending on the number of observation epochs for that 
star. 

2.2. Initial Sample 







































































































Runaway M Dwarf Candidates from SDSS 


3 


Data are taken from the Wll spectroscopic catalog, 
which consists o f 70,841 visually inspect ed M dwarfs 


from the SDSS (York et al. 2000) DR7 (Abazajian et 
al. 2009). The proper motio ns m the Wll catalog are 


based on Munn et al. (2004), which omits proper mo¬ 
tions > 300 mas yr“^. We thus retrieved pr oper motions 
for 47 h igh proper motion stars using the [Theissen et| 


aL]|2015| (submitted) catalog, in which combinations of 
SDSS, WISE, and 2MASS baselines determined updated 
proper motions for these RdM candidates. 

Magnitu des in Wll are extinc tion-corrected using the 
method of Schlegel et al. (1998). These corrections rely 
on column densities calculated from 100 and 240 /rm 
Galactic emission along the line-of-sight out to infinity, 
and therefore may overestimate t he true extinction to 
nearby stars. Jones et al. (2011) provides a more ro¬ 
bust Galactic extinction map by htting extinction curves 
to the fluxes of ~56,000 M dwarfs in the spectral range 
5700-9200 A from the Wll catalog. We correc t for ex¬ 
tinction using the method of 
for M dwarfs not present in t 


J ones et al.| (|2011|), e xcept 
|Jones et al.| (|201ip cata- 
log, where w e apply the extinction correction s of Schlegel 
et al.|(|1998|) . Extinction corrections from thelJones et al. 

luni Tk' _ r cf? act _ 1 1 'L _ _ 


(|2U11 ) catalog were available for 57,067 stars! 'i'he mean 


extinction correction for our final sample is 0.18 magni¬ 
tudes. Given that the mean apparent magnitude of our 
sample is (r) = 17.1, and the mean distance from the Sun 
is 686 pc, our stellar distances would have been offset by 
^8.6% if extinction corrections using either catalog are 
neglected. Most of our other uncertainties (e.g., metal- 
licity, spread in the photometric parallax relationship), 
are typically larger than this, so the effect of extinction 
on distance is minimal in comparison. 

2.3. Distances and Kinematics 

Because trigonometric parallaxes are not available for 
many faint stars, such as M dwarfs, alternate relations 
must be employed. (Bochanski et al. 2010) derived 
a photometric parallax relation between Mr and the 
extinction-corrected r—z color of the star. Given two dif¬ 
ferent hlters, the color of a star determines its distance in 
a particular band A using the distance modulus relation 


mAi - Mai(wai - rnx^) = 5 logd- 5, 


( 1 ) 


where d is the distance, m\^ is the apparent magnitude 
in one filter, mvi — nrix^ is the color from two filters, 
and Mxi is the absolute magnitude, which is a function 
of the color. Using a sample of nearby stars with known 
absolute magnitudes derived from trigonometric parallax 
relations, Bochanski et al. (2010) derived the photomet¬ 
ric parallax relation 

Mr = 5.190-b2.474(r-z)-b0.4340(r-^)^-0.08635(r-z)^ 

( 2 ) 

and the residual rms scatter in the relation is 0.394. Dis¬ 
tances in the Wll catalog were calculated using this pho¬ 
tometric parallax relati on. 


As first described by Sandage & Eggen (1959), lower- 
metallicity stars will have higher surface temperatures 
and thus higher luminosities than higher-metallicity stars 
of the same mass. Gorrections to the absolut e magnitude 
may be d erived using iron abund ance (e.g., Ivezic et al. 
2008). As Bochanski et al. (2013) demonstrate, the pho- 
tometric parallax relation must be modified for the lower 


Table 1 

M Subdwarf Model Parameters®' 


Class 

r — z 

Mr 


sdM 

1.05 

10.0996 

0.422153 


1.35 

10.5971 

0.422138 


1.65 

11.5888 

0.422124 


1.95 

12.9130 

0.422156 


2.25 

13.3980 

0.422160 

esdM 

1.05 

9.98573 

0.421217 


1.35 

11.4767 

0.420383 


1.65 

12.5346 

0.421741 

usdM 

1.05 

11.1409 

0.415195 


1.35 

12.1264 

0.413827 


Figure 4a of [Bochanski et al.| | |2013^ 


metallicity classes (subdwarfs, extreme subdwarfs, and 
ultra subdwarfs). By taking metallicity into account, 
which we do in this study, the distance uncertainties 
are actually reduced. However, precise metallicities for 
many M subdwarfs have not been measured, and instead 
proxies for metallicities have been developed to estab¬ 
lish metallicity classes. These classes are parameterized 
by the metallicity-dependent quantity C, which is a rela¬ 


tion between the GaH and TiO molecular indices (Dhital 


et al. 2012). The indices upon which C are based have 


been employed as t emperature and metall icity indicators 
for low-mass stars (Bochanski et al. 12007 and references 

lOT 


therein). The parameter Q is calibrated over spectral 
types M0- M3, but suffers fro m large spreads near solar 
metallicity (|Woolf et al.|2009|). |Lepine et al.| (|2007|) clas¬ 
sified stars as subdwarfs, extreme subdwarls, and ultra 
subdwarfs, for C, in the ranges 0.5 — 0.825, 0.2 — 0.5, and 
less than 0.2, respectively. 

For stars with ^ > 0.82 5 (dwarfs), we use the photo¬ 
metric parallax relation of Bochanski et al. (2010), which 
provides distance uncertainties of ~207o due to scatter 
along the main sequence. This is a conservative estimate, 
since no metallicity information was used to derive the 
photometric parallax relation. For stars of lower metal¬ 
licity classes, we ad opt empirical models of Mr vs. r — z 
using the results of Bochanski et al. (2013). Specihcally, 
we ad opt the models in h'lgure 4a of [Bochanski et ^ 
(2013), which presents separate empirical fits for subd¬ 
warfs, extreme sub dwarfs, and ultra subdwarfs. Table 
fT] presents the parameters that Bochanski et al. (2013) 
aerived for their M subdwarf niodeis (J. J. Bochanski, 
private communication). 

The discrepancy in Mr between metallicity subclasses 
is ~1 magnitudes, which for a subdwarf with apparent 
magnitude r = 20, translates to a distance d iscrepancy 
of ^0.2 kpc. Bochanski et al. (2010| 2013) find that 
the scatter in their parallax relations is (Jm ~ 0.4. To 
simulate the distance errors due that would result from 
variation in we perform Monte Carlo simulations with 
1000 realizations on r, r — z, and the extinction values, 
assuming normally-distributed errors. Our final sample 
has a mean distance uncertainty of ~26%, primarily due 

to variation in _ 

The RVs, d istances, and pr oper motions fro m |Munn| 


et al. (2004 or if missing, Theissen et al. submit¬ 
ted) were used to calculate the T, F, W velocities of each 
star, where U points from the Sun toward the Galac¬ 
tic center parallel to the midplane, V points in the di¬ 
rection projected from the Sun’s velocity onto the mid- 
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lowing the convention of Zheng et al. 

(2014) 

we adopt 

250 km s ^ for the velocitv of the Local Standard of Rest 

(LSR; Reid et al. 2009 McMillan & Binney 

2010) and 

(C/o, Vb, Wo) = (11.1,12.24,7.25) km s"" for t 
motion of the Sun with respect to the LSR 

le peculiar 

(Schonrich 


subdwarts, extreme subdwarfs, and ultra subdwarfs in 
our definition of runaway M dwarf candidates (RdMs). 

For the purpose of constructing a high quality sample 
of RdMs from Wll, we kept stars that had flags set to 
GOODPHOT = 1, GOODPM = 1, and WDM = 0. We also relaxed 
the criterion on nFit and include stars down to nFit 
= 4. We excluded stars whose S/N of the spectrum near 
Ha was less than 3. To ensure that our sample contains 
stars with low proper motion uncertainties, we also culled 
from our sample all stars farther than 1 kpc from the 
Sun. Based on the provided distances, proper motions, 
and RVs, we constructed a reduced sample of 196 RdMs 
whose GC velocities exceed 400 km s“^. 

2.4. Spectra 

To confirm the RVs for the RdM candidates, we in¬ 
spected two sets of sodium absorption transitions present 
in the spectra of the RdMs. Appendix presents the 
sodium transitions under consideration, and their rest 
frame wavelengths in vacuum. We visually examined 
the Doppler shift of these sodium transitions to validate 
the RV measurements. Heliocentric RVs and their un¬ 
certainties, recorded in Wll, were determined using the 
method of B07b. Some of our RdMs had RV uncertain¬ 
ties in the Wll catalog exceeding 100 km s“^. We thus 
applied an alternative cross-correlation function (CCF) 
routine to estimate an additional RV uncertainty. We 
first determined the spectral type (e.g., MO) of each RdM 
and splined its corresponding template to the SDSS res¬ 
olution. We cross-correlated the SDSS spectrum with 
the template in five wavelength windows and splined the 
cross-correlation function (CCF) to determine the lag at 
which the CCF is a maximum. The windows range from 
6531 to 7760 A in step sizes of 0.0150 log (A), or ~250 A. 
Our justification for these windows is that M dwarf SDSS 
optical s pectra tend to hav e higher S/N ratios toward the 
red end (Jones et al. ]2MT| ), which produces smaller vari¬ 
ations in the lag. We calculate the standard error of the 
lags about their mean to determine the RV uncertainty. 

We retained stars with clearly discernible sodium tran¬ 
sitions in their spectra and RV discrepancies within a 
2(7 RV uncertainty window. As a result of our vi¬ 
sual inspection, we arrived at a refined sample of 33 
RdMs. Figure [J shows an example of the SDSS op¬ 
tical spectrum ot J184314.08-1-412258.9, an MO, whose 
GC velocity is 427.2 ± 107.8 km s“^ and whose RV is 
—317.1 ± 6.6 km s“^. Figure also presents a close-up 
of the spectrum in the vicinity of the 3p—>3d Na triplet 
(8185.51, 8197.05, and 8197.08 A). 

2.5. Proper Motions 

To investigate the robustness of our proper motions, 
we compared the proper motions for ou r RdMs from 
SDSS -FUSNO- B (used in Wll), PP MXL (|Roeser et al. 
2010|, LSPM (iLepine & Sharal20051), and a master cat- 
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Figure 1. Close-up of the SDSS optical spectrum of 
J184314.08-1-412258.9, an MO, with Galactocentric velocity 427.2di 
107.8 km and radial velocity —317.1 di 6.6 km s“ , in the 
vicinity of the 3p^3d Na triplet. The dotted gray lines located 
at 8185.51, 8197.05, and 8197.08 A indicate the rest frame tran¬ 
sitions in vacuum. The solid blue lines and the neighboring solid 
black lines indicate the expected Doppler-shifted wavelengths of 
the transitions and 2cr uncertainties in the RV, respectively. 

submitted) , which combines dat a from WISE dWrigh t 
et al. 2010), SDSS, and 2MASS (Skrutskie et al. 2006p . 


For LSFM proper motion errors, we adopt the quoted 
velocity dispersio n based on SUPERBLI NK astrometry 
(see Section 5.2 of Lepine k, Shara ( 2005| )) to give 2cr er¬ 
rors of {fia = 14.6, fis = 12.6) mas yr“^. In making our 
comparisons, we considered whether the SDSS-I-USNO-B 
proper motions shared a 2a overlap with the proper mo¬ 
tions of the other catalogs. The 1T/S'i?-|-SDSS-I-2MASS 
baselines are tied to the ICRS, whereas five of the six 
SDSS-FUSNO-B epochs use relative positions and proper 
motions. From our sample, we removed four RdMs whose 
SDSS-FUSNO-B proper motions were inconsistent with 
those of IF/S'i?-|-SDSS-|-2MASS. LSPM proper motions 
were available for 13 of the remaining 20 RdMs. Table 
1^ presents the measured proper motions and 2cr errors 
based on the various catalogs under consideration. Table 
shows that the proper motions are generally consistent 
among the catalogs to within the uncertainties. 

3. RESULTS 
3.1. Properties 

Our final sample consists of 20 RdMs with high-quality 
SDSS spectra. We present the properties of our RdMs in 
Table[^ where M is the spectral type, C is th e me tallicity- 
dependent parameter discussed in Section |2.3[ r is the 


method of 

Jones et al. (2011 or 

Schlegel et al. 1998 

if the 

KdM was not in the Jones et al. 

(21)111 catalog), V is the 

distance, KV is the heliocentric radial velocity, D, V, IT 


alog (hereafter M/lS'A-|-SDSS-F2iVIASS; Theissen et al. 


are the velocity components in the Galactic coordinate 
system, computed u s ing th e SDSS-I-USNO-B proper mo¬ 
tions of Munn et al. (2004), TV is the tangential velocity 
(proper motion scaled to distance), Utot is the velocity of 
the star with respect to the Galactic center, and fi is the 
SDSS-I-USNO-B proper motion i n right ascension (a , al¬ 
ready corrected for (leclination injMunn et al.|2004| and 
declination (S). Uncertainties in 'i’tot were calculated by 
adding the uncertainties in U, V, W in quadrature. Re- 
po rted RV uncertainties ar e the lesser of those calculated 
by Bochanski et al. (2010) or our CCF routine, outlined 
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Figure 2. Kinematic distribution of our runaway M dwarf can¬ 
didates, with la error bars in each velocity component. Our 
sample generally consists of stars with \U\ > 200 km s“^ and 
V < 0 km s“^, but no preferred speed in W. 

in Section [2^ The uncertainty in ( for RdM 3 is large; 
generally, the [Fe/H]-C relation has conside rable scatter 
at th e high-metaihcity end of the relation (Woolf et al. 
2009|). As Figure!^ shows, the RdM 5 spectrum has com¬ 
parable S/N to me spectra of our other candidates, so 
we attribute the high C. uncertainty to scatter at the high 
metallicity end of the [Fe/H]-C relation. 

Figure shows the k vs. U and k vs. W kinemat¬ 
ics of our sample, with Icr error bars in each velocity 
component. The majority of our RdMs have k veloci¬ 
ties that are moving in the opposite direction of Galactic 
rotation. With the exception of RdM 4, our sample con¬ 
sists of RdMs with U speeds exceeding 200 km s“^. Our 
sample generally consists of stars with \U\ > 200 km s“^ 
and k < 0 km s“^, but no preferred speed in W. Our 
W velocity dispersion is much lower than our U or V 
velocity dispersion. We discus s th e implications on our 
tangential velocities in Section |4.2[ 

Figure [^presents the optical spectra of our RdM sam¬ 
ple in the range 5000 — 9000 A. All of our spectra show 


clear TiO features. The extreme subdwarf RdMs 4, 8, 
17, and 19 show a significant reduction in the TiO band, 
as is expected in the spectrum of extreme subdwarfs 
(Savcheva et al. 20141. The 3s—^Sp Na doublet is par- 


ticularly prominent ih the spectrum of each RdM, while 
the 3p—>'3d Na triplet stands out less in the spectrum of 
most of our RdMs. 

3.2. Are the RdMs Bound to the Galaxy? 

To test whether each of the 20 RdMs are bound to the 
Milky Way, w e adopt the three-co mponent bulge-disk- 


halo model of Kenyon et al. (|2008 ), as it is well suited 


to study the kinematics of hypervelocity candidates. For 
GO distance r, GG radius R (that is, r projected onto 
the Galactic plane), and height 2 above the plane, the 
Galactic potential is given by 


$ = -b $ 


hi 


( 3 ) 


where the potentials for the bulge, disk, and halo are 
respectively given by 

= -GMbI (r + Ob), 

$d(i?, z) = -GMd/\J i ?2 + [od + (z2 + 52)(l/2)] 2 ^ 

^hir) = -GMh\n{l + r/rs)/r. (4) 


We adopt the parameters used by Kenyon et al. 
and update the values of the disk mass Md ana radial 
scale length Ts to be in closer agreement wit h the ob- 


( [2008t 


served m ass measurements of the Milky Way (Brown et 
al. 112014 1. The parameters of our model are given in 'i'a- 

For each star, we perform a semi-Euler (symplectic) 
integration using the 3D positions and velocities of the 
RdMs as our initial conditions, and time steps of 10® yr. 
Symplectic integration has been used to simulate the 
interaction between a binary system and the GMBH 
dBromley et al.|2006 |) for the purpose of studying hyper- 
velocity candidates (Kenyon et al.||2008D, and to model 
the orbital d ynamics of the giant p lanets in the early 
solar system (Morbidelli et al. 2007). We classify each 
RdM as bound or unbound by integrating until: 

1. the star reaches the Galact ic virial radius (Ry w) of 
250 kpc (see Section 4.2 of Brown et al.||2014 ). 


2. the star, once gaining distance from the Galactic 
center, reaches a turnaround distance within R^w, 


or 


3. a total of 2.5 Gyr (25,000 steps) elapses without 
either of the previous conditions being met. 

All of our RdMs ultimately meet either the first or the 
second condition. Table presents the results of our in¬ 
tegrations. We consider an RdM “unbound” if it reaches 
i?vir within 2.5 Gyr. For each unbound RdM, we re¬ 
port the escape time tesc, which is how long the RdM 
is expected to take to reach i?vir- Alternatively, for each 
bound RdM, we report the maximum distance i?max from 
the Galactic center. The column “% Unbound” repre¬ 
sents the percen tage of our simulations, which we de¬ 
scribe in Section 13.31 in which the RdM reaches i?max- 
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5000 6000 7000 8000 9000 

Wavelength (A) 


Figure 3. Spectra of our runaway M dwarf candidates. The spectra have been normalized and multiplied by a constant. Another constant 
has been added to each plot to stack spectra accordingly. Our RdMs show clear TiO and Na features, particularl y the 3s^3p Na doublet, 
in their spe ctra. The feature at ~5600 A is due to the dichroic cutoff of the red and blue cameras and is not real (Sako et al.|2005 Morgan 
|et al.|2012[ |. -- - 
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Figure 3. Continued. 


Table 4 

Adopted Parameters of the Galactic Potential 



dh 

Md 

a-d 

bd 

Mh 

rs 

(M©) 

(kpc) 

(M©) 

(kpc) 

(kpc) 

(M©) 

(kpc) 

3.76 

0.1 

6 X 10i° 

2.75 

0.3 

1012 

20 


3.3. Simulations 

Table shows that at least two RdMs are unbound, 
as at their current positions and velocities, they will 
reach the Galactic virial radius within 1 Gyr. To ex¬ 
amine the extent that the position and velocity affect 


the candidacy of each RdM as bound or unbound, we 
performed Monte Garlo simulations with 1000 realiza¬ 
tions for each of the RdMs. The results of our sim¬ 
ulations depend on the velocity error distribution, of 
which proper motion errors contribute most significantly. 


The p roper motion error distribution of Dong et al. 


(20111, which is based on observations of quasars, is 
modeled by a Gaussian co re plus a wing func tion. Ac¬ 
cording to Figure 3 from |Dong et al.| (|2011|), the dis¬ 
crepancy between the normalized core-l-wing distribution 
and a normalized core (exclnsively) for proper motions 
within 10 mas yr“^ decreases for fainter sources down to 
g < 20.5. Since our RdMs have Icr proper motions errors 
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Table 5 

Galactic Potential Test Results for our Runaway M Dwarf Candidates 


RdM 

Name 

U 

(km s“^) 

u 

(km s~^) 

W 

(km s“^) 

Reaches 
R ■ 7 

-ttvir • 

^esc Ol' [Rmax]^ 
(Myr or [kpc]) 

Unbound'^ 

(%) 

1 

J003012.87-184446.9 

646.0 ± 143.0 

-374.9 ± 185.3 

-24.8 ±36.5 

yes 

628.4 

77.2 

2 

J023510.64-I-004924.9 

440.5 ± 115.2 

-616.9 ± 179.3 

-158.5 ± 112.7 

yes 

854.0 

65.0 

3 

J035310.50-004928.1 

201.4 ± 121.0 

-603.4 ± 183.3 

-100.2 ± 120.5 

no 


44.0 


23.4 

4 

J100557.07-I-344549.0 

-370.9 ± 126.5 

-430.8 ± 113.8 

-208.1 ± 94.5 

no 


63.5 


24.4 

5 

J110252.68-I-274203.7 

337.4 ± 73.7 

-484.6 ± 176.1 

-35.4 ± 19.3 

no 


40.5 


11.0 

6 

J111825.99-I-090229.4 

-422.4 ± 154.7 

-466.0 ±83.9 

-110.9 ±65.9 

no 


80.1 


32.3 

7 

J120525.80-I-403509.8 

-373.2 ± 124.6 

-469.3 ±88.5 

41.8 ±29.5 

no 


49.9 


14.8 

8 

J121441.21-1-414924.8 

-254.8 ± 164.3 

-723.8 ± 154.3 

-45.5 ±44.8 

no 

[146.6] 

43.6 

9 

J131702.01-I-382435.2 

-533.7 ± 122.2 

-350.7 ± 88.8 

7.2 ± 19.9 

no 

[156.9] 

43.6 

10 

J140921.10-1-370542.6 

-238.2 ± 178.9 

-644.7 ± 165.4 

130.1 ± 69.6 

no 

[72.0] 

36.4 

11 

J142546.68-I-082717.2 

338.0 ± 73.7 

-629.2 ±97.4 

-84.5 ± 44.0 

no 

[100.7] 

28.0 

12 

J153737.72-005608.7 

-106.0 ± 19.0 

149.5 ±28.9 

96.3 ± 19.0 

no 

[39.4] 

0.0 

13 

J182547.07-I-641355.5 

-533.9 ± 153.8 

-159.9 ±30.1 

-182.8 ±48.8 

yes 

1368.9 

49.8 

14 

J184314.08-I-412258.9 

-393.9 ± 91.3 

-257.5 ±27.7 

165.1 ± 50.0 

no 

[44.8] 

6.7 

15 

J203803.39-I-145300.2 

485.6 ± 76.9 

-259.4 ± 60.9 

-70.7 ± 74.3 

no 

[69.9] 

13.7 

16 

J221152.95-I-002207.7 

431.0 ± 84.4 

-386.0 ± 103.9 

-48.4 ± 73.5 

no 

[53.2] 

9.8 

17 

J224403.71-1-231532.4 

-487.0 ± 91.3 

-314.0 ±28.6 

-243.1 ± 51.3 

no 

[173.6] 

41.3 

18 

J231405.61-I-230120.2 

-434.1 ± 124.5 

-304.1 ± 49.3 

-27.1 ± 54.3 

no 

[50.3] 

14.3 

19 

J232541.30-I-000419.6 

-366.6 ±61.1 

-312.8 ±29.1 

-147.9 ±24.7 

no 

[36.7] 

0.2 

20 

J235459.63-004133.2 

371.4 ± 16.6 

-60.8 ± 13.6 

20.0 ±9.1 

no 

[40.4] 

0.0 


^ The escape time (tesc) is reported for unbound RdMs, whereas the maximum distance (Rmax) is reported for bound RdMs. See 
Section El for further information. 

^ “% Unbound” represents the percentage of our simulations (described in Section|3.3[> in which the RdM reaches Rmax- 


within ~3 mas yr“^, we have simply adopted a Gaussian 
distribution of the SDSS+USNO-B proper motion un¬ 
certainties for our analysis. The column in Table [^titled 
“Unbound” represents the percentage of our simulations 
in which the RdM reaches i?max- 
Table E shows that RdMs with GC velocities under 
570 km s“^, are bound in >50% of our simulations. 
All but three bound RdMs are expected to turn around 
within 100 kpc of the Galactic center. About half of our 
simulations show RdM candidate 13, with GC velocity 
571.5± 164.1 km s“^, as being bound to the Galaxy. For 
our hypervelocity candidates, Candidate 1, with GC ve¬ 
locity 658.5 ± 236.9 km s“^, is unbound in 77% of our 
simulations. We suggest that this RdM may be a true 
hypervelocity star with over Itr confidence. Candidate 
2 may also be potentially unbound, although we cannot 
claim much more than ^Icr confidence due to our high 
uncertainties, of which proper motion and metallicity un¬ 
certainty contribute the most significantly. 


4. DISCUSSION 
4.1. Selection Effects 

We have presented a sample of 20 RdMs which consists 
primarily of bound tangential velocity outliers and three 
likely h yper velocity candidates above Icr confidence. In 
Section |2.3| we outlined a process by which we created 
a reduced sample of 192 candidate RdMs prior to ex¬ 
amination of their spectra and proper motions. Section 

S describes the second stage of the reduction process, in 
- ch w e examined the spectra of the RdMs by eye. Sec¬ 
tion ^ describes the third and final stage of the reduc¬ 
tion process in which we compared the proper motions 
of the RdMs among multiple catalogs. 

Table [^summarizes the number of dwarfs in each spec¬ 
tral type at each stage of the reduction process. Most 
of the later-type RdMs were removed from our sam¬ 
ple after the first stage as they tended to have noisy 
spectra and/or large RV errors. We coincidently ended 


Table 6 

Number of Dwarfs in Each Spectral 
Type at Various Reduction Stages 


Stage 1“ Stage 2^ Stage 3'^ 

MO 

62 

16 

16 

Ml 

50 

5 

4 

M2 

18 

0 

0 

M3 

21 

0 

0 

M4 

19 

0 

0 

M5 

11 

2 

0 

M6 

6 

1 

0 

M7 

5 

0 

0 


^ 192 RdM candidates, Sectio n |2^ 
^ 24 RdM candidates, Sectionji^ 

^ 20 RdM candidates, Section|2.5| 


up omitting the M5 and M6 candidates due to mis¬ 
matched proper motions between SDSS-fUSNO-B and 
1F/S'i?-|-SDSS-|-2MASS. Our selection process does not 
necessarily mean that all later-type (>M4) dMs have bad 
spectra. Out of the 2 x 10^ later-type dMs in Wll, 
c.^1600 have GC velocities exceeding 400 km s“^. Ini¬ 
tial U, V, W veloci t ies co me from proper motions taken 
from Munn et al. (20041, which are also biased toward 
ear her-type dMs. 

To test whether the process of flagging stars with bad 
spectra or Na shifts biased the kinematics of our sam¬ 
ple, we performed independent one-way ANOVA tests for 
each of U, V, W between our initial sample of 192 RdMs 
and our final sample of 20 RdMs. Table [^presents the re¬ 
sults of our ANOVA tests, where SD represents the stan¬ 
dard deviation, df is the number of residual degrees of 
freedom, F is the F-ratio, and p is the p-value. We mark 
statistically significant p-values < 0.05 with an asterisk. 
Table shows that there is no statistically significant dif¬ 
ference in the U or W kinematics between the two sam¬ 
ples; however, the final sample has a significantly lower 
(V) than the initial sample ((U) = —385.0 km s“^ com- 
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Table 7 

ANOVA Test of Sample Kinematics 



n 

M (km s 

SD (km s-i) 

df 

F 

P 

u 

192 

-49.8 

486 

210 

0.01 

0.91 


20 

-63.2 

412.5 




V 

192 

-57.1 

444.2 

210 

10.61 

< 0.05* 


20 

-385.0 

212.8 




w 

192 

203.9 

703.4 

210 

2.62 

0.11 


20 

-51.4 

109.4 





Statistically significant p-values are marked with an asterisk. 


pared to {V) = —57.1 km s ^). A similarly negative {V) 
for subdw arf RdMs has been p reviously observed (see 
Table 7 of Savcheva et al.|[2014| {V) = —458.3 km s“^). 
Negative V velocities are consistent with asymmetric 
drift and dynamical interactions. We do not see any 
obvious reason for why our selection process, based on 
inspection of their spectra by eye and comparison there¬ 
after of their proper motions, would significantly offset 
the V kinematics of our final sample. 

We thus ask whether or not we should expect to ob¬ 
serve an isotropic M dwarf stellar velocity distribution 
near the Solar neighborhood. If so, then the TV should 
be \/2 times larger than the RV. Figure|^shows that our 
sample of 20 RdMs ten ds toward higher TVs than RVs. 


Palladino et al. (2014a I noted that their sample of 20 G 
and K dwarfs suhered fr om a high TV to RV ra tio (see 
Figure 1 in the erratum Palladino et al. 2014b), which 
they suggest is the result of possible contamination due 
to large proper motion errors. Their Monte Carlo analy¬ 
sis adopts th e proper motion error distribution of |Dong| 


et al. (2011), which models the proper motion errors of 
quasars by a Gaussian core plus a wing function. The 
analysis by Palladino et al. (2014a) revealed that 13 of 
their 20 G and K dwarfs retain hypervelocity status, and 
all of their candidates have less than a 25% probability 
of being interlopers. In Table we checked that the 
proper motions of our sample are generally consistent 
among multiple catalogs to within the uncertainties, so 
we suggest that the likelihood of contamination due to 
proper motion errors in our sample is minimal. 

Since the SDSS footprint is largely out of the Galactic 
plane, we expect the W velocity should trace the RVs. 
Our mean W and RV sample velocities are —54.8 km s“^ 
and —65.1 km s“^, respectively, which thus sug gests an 
overal l low R V contribution to the total velocity. Kenyon] 


et al. (2014) compared their simulated distance-limited 
proper motion and RV distributions for hypervelocity 
stars and runaway stars in the range 1 — 3 M(ti for both 
su pernov a and dynamical ejection scenarios. jKenyon etj 
al (2014|) found that while the RVs initially published by 


Palladino et al. (2014a) are consistent with their model 
predictions, the large proper motions are not. Using im¬ 
proved proper motion measurements for 14 of the 20 
stars based on all available images i ncluding SDSS at 
different epochs, Ziegerer et al. (2015) showed that 11 of 
these 14 stars are bound in over 997o of their simulations. 


Ziegerer et al. (2015) concluded that 13 stars belong to 
the disk population and have kinematics consistent with 
a bin ary supernova or pe rhaps a disk runaway mecha¬ 
nism (Kenyon et al. ]M4l ), while the remaining star may 
belong to th e halo. 'Fhe models developed by jKenyonj 


et al. (2014), however, have a lower mass limit of 1 Mq, 


o 



Figure 4. Tangential vs. radial velocities of our RdM sample, 
with la error bars. The red lines correspond to a perfectly isotropic 
velocity distribution, where the tangential velocity is \/2 times 
larger tha n the radial velocity. As in Figure 1 of [Palladino et| 
jal.jl 2014b|l, the majority of our sample has much higher tangential 
velocities than radial velocities. Our proper motion comparisons, 
in Table 1^ however, suggest that our sample is not contaminated 
by high proper motion errors. 

and so their predicted velocity distributions may not be 
reliably applied to the lowest mass M dwarfs, ~0.1 Mq, 
which are the dominant stellar constituent. 


4.2. Tangential Velocities 

One initial observation of the stellar kinematics in Ta¬ 
ble is that the W velocity dispersion, 109.4 km s“^, 
is much lower than the U and V velocity dispersions, at 
486.0 km s“^ and 444.2 km s“^, respectively. The SDSS 


footp rint is largely out of the Galactic plane (York et al. 
20001, so we should expect a low W dispersion to be as- 


sociated with a small RV dispersion. If we compare the 
absolute values of W and RV over the sample, then we 
calculate a total W to RV ratio of 0.89. Likewise, If we 
compare the absolute values of U, V, and TV over the 
sample, then we compute a total U to TV ratio of 0.66 
and a total V to TV ratio of 0.68. 

To assess whether the tangential velocities suffer large 
systematic errors, we compare the kinematic dispersion 
with that of 14 M subdwarfs wi th GG velocities exceed¬ 
ing 525 km s“^, published in Savcheva et aj^ ( 2014 ) 
and the G an d K dwarf sample, published in jra 


( 

2014), 

■^a 

ladino 


et al. (2014a) and Ziegerer e ](_aj^ 20]_^^^__^The veloc- 


are 


ity dispersions au,av, cryy m Savcheva et al. (|2014 1 
381 km s“^, 283 km s“^, and 191 km s“ , respec- 
tively, which suggests that a low velocity dispersion in 
the W component is not u nique to our study. More¬ 


over, 


Savcheva et al. (2014) observed lower W velocity 


dispersions than U or V velocity dispersions for each M 
subdwarf metallicity class in their samp le, as well as for 


non-hypervelocity stars in their sample (Bochanski et al.j 
2013). Fuchs et al. (2009) observed a sirnilar trend tor iVl 

dwarfs in the thin disk from SDSS DR7^ _ 

For c o mparison with t he stu dy of jPalladino et al.| 
(|2014a|), jPalladino et al.j (|2014a|) used on ly SDSS pho¬ 
tometry, whereas Ziegerer et ah (2015) derived new 
proper motions for 14 ot 20 stars “by combining all 
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astrometric information at hand from digitized plates 
and modern surveys.” Their published proper motions 
give velocity dispersions au,av,crw of 59.5 km s“^, 
43.6 km s“^, and 91.6 km s“^, respectively. The low 
radial velocities of their sample are consistent with a low 
W velocity dispersion. The observed ^ and F velocity 
dispersions observed by Ziegerer et al. (2015) are lower 


than the W velocity dispersi on. 


Why does Savcheva et al. (2014) find lower W veloc- 
ity dispersions t han U or V velocity dispersions, when 
Ziegerer et al. (2015) find the opposite? The answer 
may lie m the sampled kinematics of the stars. The 
mean Galact ocentric velocity of th e 14 stars under con¬ 
sideration in Ziegerer et al. (2015) is 282 km s“^, and 
only 3 stars retain Galactocentric velocities exceeding 
400 km s~^. The 14 M subdwarfs with published in 


Savcheva et al. (2014) have a mean Galactocentric veloc¬ 
ity of 674 km s“^, which is along the tail end of stellar 
kinematic models inv olving supernova eje ctions and dy¬ 
namical interactions (|Kenyon et al.||2014l and references 
therein). Our proper motions are confirmed by other 
catalogs (Table and while our distance uncertainties 
are larger, they cannot account for the high tangential 
motions that we observe. Since the mean Galactocentric 
velocity of our RdM sample is 486 km s“^, it is perhaps 
to be expected that our velocity dispersions are not i n 
a greement with those ob s erved by Ziegerer et al. (2015). 


Sellwood and Binney| (2002) identihed a mechanism 
that explains the radial migration of stars due to ex¬ 
changes in angular momentum caused by transient spiral 
modes. In these modes, transient spiral arms can move 
corotating stars without heating the disk. These modes 
are consistent with a break in the stellar surface density 
that propag a tes radially outward w ithin the disk |Roskar| 
(2008). Vera-Giro et al. (2014) simulated the evolu¬ 
tion or a galaxy representative of the Milky Way and the 
effect of the spiral modes on the stellar velocity disper¬ 
sion. Our M dwarfs have systematically negative V ve¬ 


locitie s , so a 5 Gyr evolution, as simulated by |Vera-Ciro| 
et al. (2014), might be a long enough time for pertur¬ 
bations due to spiral arm torques to perturb the orbits 
of these stars so that they lag b e hind the rest of the 
disk. Figure 5 of|Vera-Ciro et al. (2014) show that the 
Galactocentric radial and azimuthal velocity dispersions, 
which represent au and ay at the location of the Sun, 
naturally increase ^30% more than the vertical velocity 
dispersion, aw- Tbe increase interestingly seems to peak 
particularly near the Galactocentric radius of the Sun of 
^8 kpc. The higher au, ay compared to aw is consis¬ 
tent with the kinematics of our RdM sample and of the 
hypervelocity star sample published in Savcheva et al. 
(2014). 

Dynamical heating has been alter natively suggested 


Savcheva et al. 20141). In particular, Savcheva et al. 


a source for M dwarf disk runaways ( West et al. 2008 


examined the uD velocities of field M dwarfs from 
/fl and 3517 new subdwarfs and argued that thin-disk 
dynamical heating may cause the older stars to have 
highly elliptical Galactic orbits that take them out of 
the disk. These orbits would allow for higher observed 
velocity dispersions in the solar neighborhood. Simula¬ 
tions of dynamical heating predict typical W dispersions 
of < 20 km s“^ for M dwarfs (see Figure 6 of West et al. 


2006| and Figure 5 of|West et al.||2008|). However, these 
studies do not target hypervelocity star candidates, and 
so we do not quantify whether such models may be gen¬ 
eralized to explain the high velocity dispersions in our 
sample. 

4.3. Speculation of Origins 

To examine the extent to which the RdMs share a GG 
origin, we calculate the absolute difference in the angle 
between the GG position vector in the X, Y, Z coordi¬ 
nate syste m an d the GG velocity vector, Wtot (defined 
in Section |3.1l). In Figure we plot, for each RdM, 
the absolute difference (axis label |A0|) between the GG 
position angle and Vtot, with Icr error bars calculated us¬ 
ing the procedure in Appendix The values A9 — 0° 
{A9 = 180°) represent motion directly radially out from 
(in toward) the Galactic center. The value A9 = 90° rep¬ 
resents circular motion around the Galactic center. Our 
RdMs, which are all within ~1 kpc of the Sun, tend not 
to exhibit circular motion around the Galactic center, as 
exhibited by the absence of RdMs with A9 = 90°. 

Figure [^presents a case against a GG origin for about 
half of our RdMs. The range of possible differences in 
GG position angle and GC velocity angle rules out the 
values 0° (radially outward from the Galactic center) 
and 180° (radially inward toward the Galactic center) 
to within la for 12 RdMs and 2a for 8 RdMs. Hence, 
only about half of our RdMs have kinematics consistent 
with a GG origin. We note that some RdMs may not ac¬ 
tually reach the Galactic virial radius and instead have 
turned around. Interestingly, RdMs 1 and 2, which have 
an inward velocity component, are likely unbound, as es¬ 
tablished in Table yet they are getting closer to the 
Galactic center. Sum an observed velocity is inconsistent 
with an ejection mechanism from the center of the Milky 
Way, and so the origin of their high tangential veloci¬ 
ties is unlikely to be Galactocentric. Other acceleration 
mechanisms, the origins of which we describe below, may 
influence an otherwise exclusive GG origin. 

One possible alternate ejection mech anism is through 


a supernova explosion (Tauris 2015 and references 


therein), as the energy of the ejection could push the 
companion star forward. In the s upernova eject ion sce¬ 
nario pertaining to G and K stars, Tauris (2015) showed 
that GG velocities up to ~770 and ~1280 km s“^ are 
possible for these stars, respectively. Such an explosion 
would also contribute metals to the companion star. As 
16 of our 20 RdMs have f > 0.825 (they are dwarfs, 
rather than a type of subdwarf), we cannot rule out the 
possibility that one or more of our RdMs may have been 
accelerated from a supernova. Since abundance measure¬ 
ments for dMs are uncertain to ~0.2 dex, we do not have 
any way to probe if the RdMs are being enriched given 
current abundance measurements. 

Another plausible mechanism for accelerating main- 
sequence sta rs is dynamic ejection through multi-body 
interactions. Kollmeier et al.| (2009) initially determined 

for old popu- 
U sing a Solar- 


a maximum ejection rate of 35 Myr“^ 
lation stars from the Galactic center, 
neig hborh ood survey volume of 265 kpc^, Kollmeier et 


al. (2010| found that ejection rates are even lower for 


late-type stars and stars with sub-solar metallicities. 
Our RdM sample is systematically metal-rich and so is 
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RdM 


Figure 5. Absolute difference between the GC position angle and u^ot (vertical axis, |A0|), for each RdM, with Icr error bars. RdMs with 
|A0| < 90° are moving away from the Galactic center, whereas RdMs with |A0| > 90° are moving towards the Galactic center. About half 
of our RdMs have kinematics that are not consistent with a strictly radial origin out away from, or toward, the Galactic center. RdMs 2, 
4, and 14 are expected to leave the Galaxy, but are currently getting closer to the Galactic center, which suggests that they may not have 
been ejected out of the center of the Milky Way. 

ejection rates out of M32. The predicted space density 
of M31 hypervelocity stars of all types near the Sun is 
^0.001 kpc^. Given our survey volume, the probability 
of discovering a single M31 hypervelocity star, based on 
its predicted space density near the Sun, is low (^0.1%). 
For comparison, we observed 20 RdMs within our survey 
volume, so to first order, the probability of any of our 
RdMs comi ng from M31 is 1 — (1 — 0.004)^°, or ~8%. 


at lea st consistent with the result s of iKollmeier et al.l 
(2010). iGvaramadze & Gualandris (2011) simulated the 
ettect of a three-body system consisting or two high mass 
(>80 M 0 ) stars and a flyby star in the mass range 3- 
80 M©. Ac cording to their simulations (see F igure 4 and 
Figure 9 of Gvaramadze fc Gualandris]|2011 ), about 1% 
(10%) of simulations tor each of the stars with masses 
<20 Mq had escape velocities exceeding 400 km s“^ 
(200 km s“^). 

Because these simulations include only very high mass 
stars, which are short-lived and comprise a tiny fraction 
of the stellar population, the actual probability of observ¬ 
ing such a high escape velocity due to i nteracting contact 
binaries is low. Perets & Subr. (2012) modeled the rate 
of high-velocity runaways frorri binaries in dense, com¬ 
pact clus ters or cluster cores, w ith primary star masses 
> 4 Mq. Perets & Subr. (2012) showed that the ejection 


rate of hyper-runways from binaries in dense, compact 
clusters or cluster cores is on the order of 10-20 stars 
per 100 Myr, of which only one or two such stars have 
Vesc > 400 km s“^. Our survey volume covers about one 
quarter of the sky out to 1 kpc. On the basis of these 
analyses, it is statistically unlikely that our RdMs were 
accelerated to potentially hypervelocity speeds through 
a dynamic ejection mechanism. The kinematics of our 
sample, as illustrated in Figure and Figurenonethe¬ 
less show that high TVs dominate the kinematics of the 
majority of our RdMs. 

A third plausible mechanism for generating high tan- 
gental v elocities i s thr ough accelerations from nearby 
galaxies. Li et al. (2012) proposed that some hyperveloc¬ 
ity stars may be accelerated through tidal disruptions of 
the Milky Way’s dwarf galaxies or gravitational interac¬ 
tions with the GMBH of M31 or M32. We briefly consider 
how plausible these alter native mechanis ms may explain 


the origin s of our RdMs. Lu et al. (2007) and Sherwin et 
al. (2008) respectively proposed that some hypervelocity 


In addition, Abadi et al. (2009) suggested that the dis¬ 
ruption of tidal dwarf galaxies in the Milky Way could 
have produced the stream of hypervelocity stars in the 
constellation Leo {I « 230°, b « 60°). 

We quantify whether or not our RdMs may have orig¬ 
inated from the direction of M31 or the Leo stream. We 
arrive at formulas for the uncertainty in the Galactic lat¬ 
itude b and Galactic longitude I of the velocity angle 
using the procedure in Appendix O Table [^presents the 
differences and 2a uncertainties between each RdM ve¬ 
locity direction, in Galactic coordinates, and t he direc¬ 
tion of M31 and the Leo stream. According to |SherwIn| 


stars may be ejected out of M32 and M31, with higher 


et al. (2008), for the purpose of analyzing hypervelocity 
stars, one may reasonably “assume a static configuration 
over the last 10^° yr, with M31 and the [Milky Way] 
separated by their current distance.” In Table Ai/j is 
the difference between the angle of the source (M31 or 
Leo) and the RdM velocity angle. The value Atp = 180° 
represents a star whose velocity angle is opposite the po¬ 
sition angle of the source; that is, the star came from that 
source, assuming a constant flight path. The uncertainty 
in Ai/^source is the uncertainty in the RdM velocity angle, 
ag u, which we calculate using the procedure in Appendix 

m' 

Table shows that seven RdMs have velocity angles 
that point opposite the position of M31 to within 2a. 
RdMs 1 and 2, with GG velocities of 659 km s“^ and 
595 km s~ ^, re spectively, are likely unbound, as analyzed 
in Sect ion [ 3(2 Assuming constant GC velocity and a dis¬ 
tance of 800 kpc to M31, the flight times of these stars 
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Table 8 

RdM Velocity Angles and Position Angles of M31 and Leo Stream 


RdM 

Name 

b (deg.)^ 

1 (deg.)*" 

Ai/Jmsi (deg.)*"± 

Ai/iLeo (deg.)*"’** 

1 

J003012.87-184446.9 

-2.2 ± 13.5 

349.1 ± 72.1 

127.5 ± 73.2 

106.0 ±73.2 

2 

J023510.64-I-004924.9 

-15.5 ±37.3 

320.2 ± 81.1 

138.5 ±86.7 

103.4 ± 86.7 

3 

J035310.50-004928.1 

-13.8 ± 50.6 

299.7 ± 107.4 

144.6 ± 116.4 

92.2 ± 116.4 

4 

J100557.07-I-344549.0 

-26.8 ± 7.0 

206.0 ±30.6 

76.1 ± 28.4 

89.0 ±28.4 

5 

J110252.68-I-274203.7 

-4.9 ±8.3 

325.2 ± 98.4 

144.6 ± 98.4 

96.9 ±98.4 

6 

J111825.99-I-090229.4 

-13.2 ± 13.3 

207.1 ± 2.4 

81.5 ± 13.5 

75.4 ± 13.5 

7 

J120525.80-I-403509.8 

5.5 ± 26.1 

210.4 ± 7.9 

91.3 ± 27.2 

56.5 ±27.2 

8 

J121441.21-1-414924.8 

-4.8 ± 13.4 

241.7 ± 35.8 

116.1 ± 38.2 

65.5 ±38.2 

9 

J131702.01-I-382435.2 

0.8 ±9.7 

190.7 ± 28.9 

71.3 ±30.5 

66.5 ±30.5 

10 

J140921.10-1-370542.6 

15.8 ± 70.4 

238.9 ± 44.3 

121.1 ± 80.9 

44.7 ± 80.9 

11 

J142546.68-I-082717.2 

-9.4 ± 17.3 

311.7 ± 53.2 

147.3 ± 55.2 

94.1 ± 55.2 

12 

J153737.72-005608.7 

13.1 ± 7.2 

104.9 ± 14.3 

38.2 ± 15.7 

94.8 ± 15.7 

13 

J182547.07-I-641355.5 

-18.7 ± 0.1 

170.4 ±25.1 

46.1 ± 23.8 

92.1 ± 23.8 

14 

J184314.08-I-412258.9 

22.7 ± 44.1 

181.1 ± 16.0 

73.3 ± 46.4 

50.4 ±46.4 

15 

J203803.39-I-145300.2 

-8.3 ±39.3 

358.9 ±30.0 

116.0 ± 49.1 

115.8 ±49.1 

16 

J221152.95-I-002207.7 

-6.1 ± 37.8 

342.5 ± 62.2 

130.9 ± 72.2 

106.4 ± 72.2 

17 

J224403.71-I-231532.4 

-26.3 ± 1.9 

187.5 ±8.0 

60.1 ± 7.4 

93.1 ± 7.4 

18 

J231405.61-I-230120.2 

-3.5 ± 26.3 

187.1 ± 19.2 

66.3 ±32.6 

71.8 ±32.6 

19 

J232541.30-I-000419.6 

-21.7 ± 0.6 

189.7 ± 11.7 

63.1 ± 10.9 

88.0 ± 10.9 

20 

J235459.63-004133.2 

2.7 ± 4.4 

27.0 ± 2.5 

94.9 ±5.1 

114.7 ±5.1 


^ Uncertainties are 2cr. 

M31 is located at b = -21.573°, I = 121.174°. 

^ The Leo stream is located at approximately b — +60°, I — 230°. 


are approximately 1.2 Gyr and 1.3 Gyr, respectively. Ta¬ 
ble]^ also shows that four RdMs have velocity angles that 
point along the direction of the Leo stream to within 2cr, 
which suggests that no more than four RdMs could have 
originated from the Leo stream. 

5. SUMMARY 

We have presented a sample of 20 runaway M dwarf 
candidates (RdMs) with Galactocentric (GG) velocities 
exceeding 400 km s“^ within 1 kpc of the Sun assessed 
their candidacy for being unb ound from the Mi lky Way. 


Our sample is taken from the West et al. (2011 ) catalog. 


We correc ted for extinction using the method o: 


Jones et 


al. (2011), or, i f the star was n o t pres ent in that catalog, 


the method of Schlegel et al. (1998). We updated the 


catalog distances by corre cting tor the effect of metal- 


licity using the method of Bochanski et al. (2013). Our 
RdMs pass a series of photometric and proper motion 
processing flags to minimize contamination. To ensure 
that our sample contained high-quality SDSS spectra, 
we inspected the optical spectra of each RdM by eye and 
flagged stars with noisy spectra or sodium shifts incon¬ 
sistent with their measured radial velocities. 

To ensure that our final sample was not contaminated 
by high proper motion errors, we compared the proper 
motions of the remaining RdMs a mong multiple cat alogs. 


MXL (|Roeser et al. 

20101), the Lepine snara Proper Mo- 

tion catalog (LSFivI 

, Lepine & Shara 

2005 

), and a mas- 

ter catalog combining data trom the 

Wide-tield Infrared 


Survey Explorer (WISE, Wright et al. 2010), SDSS, and 
the Two-Micron All-Sky Survey ( 2MAS S ), the analysis 


of which is being conducted separately (Theissen et al. 
submitted). We fl agged stars whose proper motions m 


West et al. (2011) were inconsistent with those of the 
other catalog^ I n total, we retaine d 20 stars from the 
original catalog of West et al. (2011) 

To determine if each RdM is unbound, we modeled 
the Galactic potential using a bulge-disk-halo profile 


(Kenyon et al. 

2008 Brown et al.||2014 

). We found that 

hve of our RdiV 

Is are expected to reach the Galactic virial 


radius. Our fastest RdM, with Galactocentric velocity 
658.5 ± 236.9 km s“^, is a possible hypervelocity candi¬ 
date, as it is unbound in 77% of our simulations. The 
majority of our RdMs are likely disk runaways or halo 
objects. Hence, our study thus provides some confidence 
in the possible discovery of one hypervelocity M dwarf. 
About half of our RdMs have kinematics consistent with 
a GC origin. Seven of our RdMs have kinematics consis¬ 
tent with an ejection scenario from M31 or M32 to within 
2 cr, although our distance-limited survey makes such a 
realization unlikely. No more than four of our RdMs 
may have originated from the Leo stream. We instead 
propose that a series of multi-body interactions within 
the Galactic disk and/or accelerations from supernovae 
are plausible alternative mechanisms. If Gaia^ meets 
its predicted astrometric performance, which would yield 
proper motion uncertainties of 0.035-0.17 mas yr“^, then 
it should unambiguously determine whether or not our 
fastest RdM is indeed a hypervelocity star. The growing 
number of low-mass tang ential velocity outlier s, such as 


those in this study and in Ziegerer et al. (2015), will also 
be of great interest in the Gaia era for investigating their 
possible origins. 
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APPENDIX 

A. SODIUM TRANSITIONS 


Table 9 

Rest Frame Sodium Transitions in Vacuum 


Transition 

^vac (A) 

3s^3p doublet 

5891.58 

5897.56 

3p^3d triplet 

8185.51 

8197.05 

8197.08 


B. UNCERTAINTY IN THE ANGLE BETWEEN THE GALACTOCENTRIC POSITION AND GALACTOCENTRIC VELOGITY 

Consider a star located at GC position r moving with Galactic rest frame velocity v measured at the location of the 
star. The uncertainty in r depends on ax, o’y, <Jz, which are the uncertainties in X, Y, Z, respectively. Let and 
r_ represent the maximum and minimum GC positions, respectively, so that 

r± = {X±ax)X + {Y±aY)Y + {Z±az)Z. (Bl) 

The angle between ry and r_ is determined by 

ry • r_ 

cos{ae,r) = ]—n-i’ (^2) 

|r+||r_| 

where ag^r is the uncertainty in the direction of r. 

Likewise, the uncertainty in v depends on the uncertainties in each component, where (dX/dt) = —U, {dY/dt) = 
Y + ulsr 7 and (dZ/dt) = W. Let and r_ represent the maximum and minimum velocities, respectively, so that 

u± = —{U ± au)X -I- (U -I- ulsr ± ov)Y + {W ± aw)Z. (B3) 

The angle between Vy and V- is determined by 

, . Vy ■ V- 

cos{ag^^) = - -r.-r, B4 

|u+||u_| 

where ag^^ is the uncertainty in the angle of ti. In the limit of small angle approximations, the total uncertainty in 
the angle between r and v, which we denote ag, is determined by adding the position and velocity angle uncertainties 
in quadrature: 


0-0 = \J {(Jg,rY + {^0,vY- 


(B5) 
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C. UNCERTAINTY IN THE GALACTIC LATITUDE AND GALACTIC LONGITUDE OF THE VELOCITY ANGLE 


Consider a star moving with Galactic rest frame velocity v. To determine the uncertainty in the Galactic latitude of 
V, we project v onto a plane that measures its speed parallel to the plane and its velocity perpendicular to the plane. 
This velocity, which we denote tif,, is given by 


Vb = + (c + ulsr)2 P + WZ, (Cl) 

where p points parallel to the Galactic midplane. The uncertainty in Vb follows a similar procedure to that in Appendix 
[B| Let Vb+ and Vb- represent the maximum and minimum projected velocities, respectively, so that 

Vb± = ViU ± auY + {V + ulsr ± cTv)^ p+iW± aw)Z. (C2) 

The angle between Vb+ and Vb~ is determined by 


cos(crb) 


Vb+ ■ Vb-^ 

kb+lkb-l’ 


(C3) 


where Ub is the uncertainty in the Galactic latitude of v. 

To determine the uncertainty in the Galactic longitude of v, we project v directly onto the Galactic plane. This 
velocity, which we denote Vi, is given by 

vi = -UX + {V + vi^sk)Y. (C4) 


Let vi+ and Vb- represent the maximum and minimum projected velocities, respectively, so that 

Vi± = -{U ± (Tu)X + {V + ulsr ± crv)y. 


(C5) 


The angle between vi^ and vi- is determined by 


cos(ctj) = 


vi+ ■ vi- 


where ai is the uncertainty in the Galactic longitude of v. 


(C6) 
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